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We study approximations of billiard systems by lattice graphs. It is demonstrated that under 
natural assumptions the graph wavefunctions approximate solutions of the Schrodinger equation 
with energy rescaled by the billiard dimension. As an example, we analyze a Sinai billiard with 
attached leads. The results illustrate emergence of global structures in large quantum graphs and 
offer interesting comparisons with patterns observed in complex networks of a different nature. 



I. INTRODUCTION 

The notion of a quantum graph is known for more 
then half a century however, an intense investiga- 
tion of these structures started less than two decades ago 
0,13 as a response to progress in fabrication technologies 
which allowed to prepare microscopic graph-like struc- 
tures. Nowadays, there is an extensive literature devoted 
to the subject; for recent reviews see 0,13 and also (3- 

An attention to quantum graphs come from the fact 
that motion on their edges is easy to describe, and at 
the same time the graph structure leads to a nontrivial 
behavior. It was shown, in particular, that even a graph 
with a small number of vertices is capable of developing 
an internal dynamics rich enough to display universality 
features that are typical for the wave-chaotic behavior 
0,11,13' It is not only a theory, the results can be checked 
experimentally in a microwave graph model IToj . 

On the contrary, properties of nontrivial large-scale 
graphs have been regarded as less interesting due to the 
expected localization of the corresponding wavefunctions. 
An indication that this belief is wrong may be seen from 
the fact that complex graph-like structures, such as sys- 
tems of interconnected neurons, display surprising pat- 
terns observed, for instance, in the visual cortex of mam- 
mals ^|. It was shown in 12] that these patterns can 
be understood as a manifestation of a Gaussian random 
field, and are in this sense analogous to patterns emerg- 
ing in two dimensional quantum chaotic systems, for in- 
stance, nodal domains in a chaotic quantum billiard [13. 

A thorough investigation of such structures on graphs 
is by no means easy. To follow the mentioned example, a 
nodal domain is a connected component of the maximal 
induced subgraph of a graph F on which a function does 
not change sign; it relates the pattern formation on F to 
nontrivial algebraic questions about graph partition, etc. 
This is probably the reason why only few mathematical 
results of this type are available at present - cf. JLl, d5j. 

In this paper we are going to show that extended 
graphs support structures similar to those known from 
two-dimensional wave chaos and Gaussian random-field 
models. We will show that the structures do not depend 
on the local set up of the graph. They reflect the in- 



fluence of the graph boundary on the wave propagation 
and interference on the graph. The flnal patterns are 
"global" in the sense that they extend over hundreds of 
graph vertices. 

Our approach is based on graph embedding into a 
Euclidean space and a convergence argument; we will 
demonstrate that wavefunctions on the graph can ap- 
proximate solutions of the respective "continuous" bil- 
liard problem. The embedding assumption is naturally 
a nontrivial restriction because not every graph can be 
regarded as a subset of a Euclidean space from which it 
inherits its metric. It applies, however, to wide enough 
class of systems and allows us at the same time to cir- 
cumvent difficulties of a pure algebraic treatment. 

The technical tool to derive the approximation result 
is a graph duality adopted from |l6j| . To make the paper 
self-contained we review this theory in the next section in 
a simplified form suitable for the present purpose. Then 
we will show how solutions to the Schrodinger equation 
can be approximated by those of a Schrodinger equation 
on lattice graphs with the energy properly rescaled. 

To illustrate the result we will analyze an example of a 
lattice graph which approximates a Sinai billiard. Since 
we want to go beyond the nodal structure and to ana- 
lyze also the phase behavior of the wavefunctions we will 
study such a also from the transport point of view, at- 
taching to it a pair of semiinfinite leads; the result will 
be compared to the "true" Sinai billiard with a pair of 
leads attached. We will compare, in particular, the prob- 
ability currents and show they are similar to each other 
provided the current on the graph is properly defined as 
a vector sum of currents at the graphs links. 



II. THEORY: A GRAPH DUALITY 

By F we denote in the following a connected graph 
consisting of at most countable families of vertices V = 
{Xj : j & 1} and edges £ = {Cjn ■ {j,n) € Ic C I x I}. 
We suppose that each pair of vertices is connected by not 
more than one link, otherwise we can simply add vertices 
to any "mifltiple" edge. The set J^{Xj) = {Xn : n € 
v{i) C / \ {]}} consists of the neighbors of Xj, i.e. the 
vertices connected with Xj by a single edge is nonempty 
by assumption. The graph boundary B consists of ver- 
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tices having a single neighbor; it may be empty. We de- 
note by /e and /j the index subsets in / corresponding 
to B and the graph interior Z V \ S, respectively. 

We suppose that F is a metric graph, i.e. that it has 
a local metric structure, every edge Cjn being isometric 
with a line segment [0,£jn]- Of course, the graph can be 
also equipped with a global metric, for instance, by iden- 
tifying it with a subset of R"^ . In general the metrics may 
not coincide, however, in the next section we will identify 
them. Using the local metric, we introduce the Hilbert 
space L'^{T) :— ®(j „)g/^ L'^{0,£jn) whose elements are 
ip = {ipjri ■ ij,n) e Ic} or simply {V'j„}; in the same 
way we define Sobolev spaces on T. Given a family of 
potentials U := {Ujn} with V,„ G L°°{0,£jn) and cou- 
pling constants a := {aj G R : j G /}, we define the 
Schrodinger operator Ha = Ha{T, U) by 

Hai^jn} { -^;'„ + Ujni^jn ■ (j, n)elc} (2.1) 

on the domain consisting of all -0 with -ipjn G H^^'^(0, Ijn) 
satisfying suitable boundary conditions at the vertices 
linking the boundary values 

i^jnij) lim V^j„(a;), — 1™, V'jn(2;) , 

^0+ a;— >0+ 

. (2.2) 

where the point a; = is identified with Xj. Specifically, 
we will work here with the so-called S coupling: at any 
Xj G V we have V'jn(j) = ^im(j) =: V'j ^ov all n,m G 
Hj), and 

E i'jnU) ^ (^ji^j ; (2.3) 

it is known that among all (non-trivial) boundary con- 
ditions which make the operator (|2.1|l self-adjoint there 
are no other with wavefunctions continuous at the ver- 
tices 's^. The particular case a = represents the most 
simple boundary conditions, called usually Kirchhoff 0, 
which we will employ in the example of Sec. lIVI however, 
for the moment it is useful to consider the more general 
situation H2.3() . Furthermore, if the boundary B ^ $ we 
assume Dirichlet boundary conditions there, 

^3=0, j&lB. (2.4) 

If F is infinite one can look not only for bound states of 
Ha but also for solutions of the equation 

Ha^ = (2.5) 

referring to the continuous spectrum. To describe the 
generalized eigenfunctions we consider in such a case the 
class D\Qc{Ha) which is the subset in V^- n)£ic -^^(^j ^jn) 
(the direct sum) consisting of the functions which sat- 
isfy all the requirements imposed at ■0 S D{Ha) except 
the global square integrability. The conditions l|2.3|l de- 
fine self-adjoint operators also if the a^-'s are formally 
put equal to infinity. We exclude this possibility, which 
corresponds to Dirichlet decoupling of the operator at 



Xj turning the vertex effectively into Nj points of the 
boundary. 

We need the decoupling, however, to state the result. 
Let Ha be the operators obtained from Ha by changing 
the conditions (|2.3() at the points of T to Dirichlet and 
denote IC :— {k : fc^ G (j{H^)}. In the particular case 
when the particle is free at graph edges, Ujn ~ 0, this set 
is given explicitly as IC :— {-nnlj^ : {j, n) G Ic, n G N+}. 
We will adopt several assumptions, namely 
(i) all the potentials of the family {Ujn} are uniformly 
bounded for (j, n) d Ic , 

(a) £q := inf{ ijn ■■ {j, n) elc}>0, 
(Hi) Lq := sup{ £jn ■■ (j, n) G /£} < oo , 
(iv) Na := max{ #i^(j) : j G / } < oo . 

To formulate the result, we need a few more notions. 
On Cnj = [0,£j„], where the right endpoint identified 
with the vertex Xj, we denote as Ujn and Vjn the solu- 
tion to — /" -I- Ujnf ^ k^f which satisfy the normalized 
Dirichlet boundary conditions 

Ujn(^jn) = l-(Uj„)'(^j„) = , Vjn{0) = l-(fj„)'(0) = ; 

their Wronskian is naturally equal to Wjn — —Vjn{£jn) — 
Mj„(0). After this preliminary we can specify the result 
of 16] to the present situation. 

Theorem: (a) Let assumptions (i)~(iv) be satisfied and 
suppose that -0 G Dioc{Ha) solves (|2.5|l for some k ^ 
ICa with fc^ G R, Irak > 0. Then the corresponding 
boundary values H2.2|l satisfy the equation 

(2.6) 

Conversely, any solution {ipj : j G /j} of the system to 
(|2.6() determines a solution of (|2.5|l by 

i^jnix) = — ^ Ujn{x) - -y- Wj„ (x) if n G l/(j) D Ij , 

if) ■ 

ijjn{x) = - —y-Vjn{x) if riGi/(j)n/B. 

Wjn 

(b) Under (i), (ii), ip G i^(F) implies that the solution 
{ipj} of the system (|2.6|l belongs to £^(/i). 

(c) The opposite implication is valid provided (iii), (iv) 
also hold, and k has a positive distance from from IC. 

III. APPROXIMATION BY LATTICE GRAPHS 

As the next step, let us inspect how the above duality 
looks under simplifying assumption: we will suppose that 

(a) all the graph edges have the same length £ > and 

(b) all the potentials Ujn vanish. Then the "elementary" 
solutions can be made explicit, 

Ujn{x) — J sinfc(x — £) , Vjn{x) = — sinkx , 
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with the Wronskian Wj, 
tern of equations (|2.6|l becomes 

'ipn — ipj cos k£ 



E 



sin ki 



^ sin k£, and the dual sys- 

)S 

+ ajVj=0, jG/; (3.1) 



Notice that this is true even if some of the iy{j) cor- 
respond to points of the boundary, because we assume 
Dirichlet condition (|2.4|l there so the corresponding ipn 's 
are zero. 

So far we worked with the local metric on F, now we 
will regard the graph as a subset in R'' and assume that 
the local metric coincides with the global one obtained 
by this embedding. We will not strive for a most general 
result and concentrate on an important particular case of 
a cubic lattice graph C = C'^{£) C M.'^ whose vertices are 
points {xj{£) = {jii, . . . ■ ji £ while the edges 
are segments connecting pairs of vertices in which values 
of a single index ji differ by one. 

Theorem: Let : M"^ ^ R be a smooth function with VV^ 
bounded. Put aj{£) :— V{xj)£ and consider the family of 
operators Ha{C'^{i),0) with ^ > 0. Suppose that for any 
fixed £ and fc £ M, the family {V'|} solves the equation 



(|3.1|) . and define a step function V'j? 



by 



■ipe{x) := V'^ if 



-£<(x- Xi), < -£. 
2 - V 2 



Suppose that the family {tpg} converges to a function 
i/j : M"^ — » C as £ — > in the sense that the quantities 

Sj{£) :— 'ipi^j) — ipii^j) behave as 



J2 ien{e)-e,i£))=o{£'). 



(3.2) 



Then the limiting function solves the equation 

-A^{x) + V{x)^{x)^iyk^i;{x). (3.3) 

Proof: Let / be a C^-smooth function, using its Taylor 
expansion to the second order we find 



f{x + £)- f{x -£)- 2 fix) cos k£ 
£fc-isin k£ 



so the right-hand side tends to f"{x) + k'^f{x) as ^ — > 0; 
in fact, the error is 0{£'^) provided f ^ . Applying 
this result to the function with respect to each of the 
V variables and combining it with the fact the family 
{i'e.ixj)} solves the equation (|3.1(l we find 

l^ipixj) + i/k'^^{xj) — V{xj)il;{xj) 
£ ^ 

-sinMj J2 ieni£)-e,i£)) + oi£) 
and the right-hand side tends to zero by assumption. 



Let us add a few comments: 

(a) The requirement k ^ K, means no restriction here, 
because for a fixed k it is satisfies if £ is small enough. 

(b) The rectangular lattice used to prove the theorem is 
not substantial. The same argument can be used, e.g., 
to prove the theorem for a graph resulting from tessel- 
lation of the plane by a lattice of equilateral triangles. 
Recall also that for rectangular lattices a similar result 
can be proven by a different method using resolvent of 
the Hamiltonian - see ^3 ■ 

(c) Notice that the limiting energy has to be rescaled to 
vk^ , where v is the dimension, roughly speaking because 
all "local" momentum components are equal. This claim 
remains valid when the we replace a rectangular graph 
with a triangular one. 

(d) The fact that the motion on the graph is locally re- 
stricted to particular directions only does not mean that 
on larger scales the particle cannot move through such 
lattice in any possible angle in a zig-zag way. To illus- 
trate the last claim recall how Fermi surface looks like on 
a 2D square lattice in the free case, aj = for any j G /. 
By IgJ it is described by the equation 

cos Oi£ + cos 6*2^ = 2 cos k£ , 

where 9i are the quasimomentum components, thus for 
small £ we have at the bottom of the spectrum 



2k^ = 9 



0{£^) 



which looks like the free "continuous" motion, apart of 
the factor 2 multiplying the energy. 

A similar result can be derived if the lattice graphs do 
not cover the whole W . Consider an open set C 
and call = Cq^{£) the subgraph of C^{£) whose vertices 
are all points Xj contained in Q.. Let 'Pq{£) denote the 
union of all closed hypercubes of C^{£), i.e. the "vol- 
ume" of such a lattice in R"^. If an edge of C^{£) belongs 
to the boundary of 'Pq{£) we delete it. It may also hap- 
pen that is non-convex, i.e. there is an axis along 
which a boundary point has neighbors in Cq{£) in both 
directions, then we regard the corresponding vertex as a 
family of disconnected vertices belonging to the boundary 
of Cq{£); we call the lattice modified in this way Cq{£). 
Mimicking the above argument, we arrive at the follow- 
ing conclusion: 

Theorem: Suppose that the potential y : ^ M is 
smooth with VV bounded and set aj{£) := V{xj)£. 
Consider the dual system associated with the family 
{Ha{C^{£),0) : ^ > 0} and its solutions {^j}. Under 
the same convergence assumption as above, the limiting 
function ^ solves the equation 



Atp(x) + V{x)tp{x) ^ vk^tpix) 



(3.4) 



with Dirichlet condition, iIj{x) = for a; G dVl. 

Let us stress an important feature, namely that the 
described result has a local character. It is especially 
important from the viewpoint of the example discussed 
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below, where we will violate regularity of the solution 
at a fixed points by attaching leads to Q. This means 
that the solution has a singularity at such a junction, 
a logarithmic one for 1^ — 2, which enters the coupling 
between the billiard and the lead. Outside the connection 
points, however, the graph approximants do still converge 
to solution of the appropriate Schrodinger equation. 



IV. EXAMPLE: SINAI BILLIARD GRAPHS 

We will consider a rectangular N x N lattice graph 
with a circular part removed reminiscent of a Sinai bil- 
liard which according to the above result such a structure 
can approximate - cf. Fig. For practical calculations 
we choose TV = 97 and aj = 0, Ujn = 0; at the graph 
boundary we impose Dirichlet conditions. The lattice 
graph spacing is set to be £ — 0.15. 




FIG. 1: Sinai billiard graph 

First we look at the nodal zone structure of one of the 
eigenf unctions - cf. Fig. 13 The vertices of the graph are 
marked as black when the value of the eigenfunction is 
positive at the vertex or white when it is negative. What 
comes out is a pattern similar to that of the hard-wall 
Schrodinger problem on the corresponding Sinai billiard. 

As we have indicated we want to compare transport 
properties of such systems in the situation when an in- 
coming and outgoing lead is attached to the graph (at the 
points (14, 40) and (59, 80)) and to the billiard at the cor- 
responding places. Adding a lead to a graph, represents 
no problems: at the incoming / outgoing vertex a semi- 
infinite is added and the resulting five edges are again 
coupled by Kirchhoff conditions, (|2.3|) with aj = 0. On 




FIG. 2: Nodal domains an eigenfunction 

the other hand, coupling a billiard to leads needs explana- 
tion. Here we use a standard method and describe the at- 
tached leads (attached antenna in the case of microwave 
billiards) by Sommerfeld radiation boundary conditions. 
Using this approach the attached lead is replaced by a 
small circle and the following boundary conditions are 
imposed on its boundary: 

^+ikiP = 2ik (4.1) 
o n 

for the incoming lead and 

^ + ik^lJ = (4.2) 
on 

for the outgoing one. The radius of the circles is much 
smaller then the length of graph bonds. We have used 
r — 0.01 for the numerical tests. Another possibility is to 
relax the regularity requirement to solution to a two- 
dimensional Helmholtz equation at the junction points. 
This approach is formally equivalent to the limit r — > 
and is described in ll9|, j20|], j2lj] and l^. However 
since we are interested in global structures that extend 
over the whole graph (billiard) the detail character of the 
connection is not important. 

Let us start with comparing the wavefunctions on the 
graph with those obtained for the corresponding two- 
dimensional billiard. A typical result is displayed on the 
Figure 131 where we have plotted the absolute value of the 
wavefunction. For the graph, on the other hand, the val- 
ues of the solution on the vertices are shown. 

Speaking of phase-related effects, a primary quantity 
of interest is the probability current which in (an open) 
billiard is given conventionally by 

j{x) = Im (^VV') (x) (4.3) 
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FIG. 3: Eigenfunction comparison in terms of probability den- 
sity. In the left picture squared modulus of the graph eigen- 
function corresponds to energy E referring to the momentum 
k = 1.65. The right picture shows the same for billiard eigen- 
function of energy 2E. The color scale ranges from zero (dark 
blue) to the maximum value (dark red). 

On the graph the current flows along the edges and 
has therefore a prescribed direction - cf. Fig. 0] - so it 
is not obvious what is the quantity to be compared to 
the two dimensional case. The natural possibility is to 
add the "horizontal" and "vertical" flows at each vertex 
by a vector summation and construct a in such a way a 
vector fleld. It turns out that leads indeed to the correct 
probability current inside the two dimensional billiard. 
We computed the currents using the above procedure are 
compared them with the current inside the two dimen- 
sional billiard that was evaluated with the help of the 
formula 14.311 . The result is plotted on the Figure [3 




FIG. 4: Probability currents on F 

In conclusion we have demonstrated the existence of 
global structures on large graphs. The structures do not 
depend on the local graph topology. We were able to 
prove that the large scale structures are the same on a 



FIG. 5: Determining the vector field. The left current is 
obtained by the vector summation while the current inside 
the billiard is plotted on the right side 



square graph, where each vertex connects 4 bond, and on 
a graph consisting of equilateral triangles when 6 bonds 
are connected at each vertex. Moreover numerical results 
show that the structures do not change also for other 
types of graphs (although a rigorous proof is missing). 

The structures extend over hundreds of graph vertices. 
They make up the manifestation of complex interfer- 
ence effects and are as such difficult to understand. The 
way out is to employ embedding of the graph into the 
appropriate ambient space and proving that the graph 
wavefunction converges to a corresponding solution of 
Schrodinger equation with the scaled energy. 
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